Recovering Entanglement by Local Operations 
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We show that any quantification of the bipartite entanglement of mixed states uniquely based 
on the density operator may lead to a paradoxical increase of entanglement under purely local 
operations. This apparent paradox is solved in the physical ensemble description of the system state 
by introducing the concept of "hidden" entanglement, which measures the amount of entanglement 
that may be recovered without the help of any non-local operation. For two noninteracting qubits 
under a low-frequency stochastic noise, we show that entanglement can be recovered by local pulses 
only. We also discuss how hidden entanglement may provide new insights about entanglement 
revivals in non-Markovian dynamics. We finally propose a simple quantum information scheme, 
implementable by all-optical setups, which gives evidence of the concept of hidden entanglement. 
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Introduction.— Entanglement, arguably the most pe- 
culiar feature of quantum mechanics, plays a key role in 
several quantum information and communication appli- 
cations, including teleportation, quantum dense coding, 
private key distribution, and reduction of communication 
complexity To work properly, all the above tasks 

generally require pure maximally entangled states. Since 
entanglement cannot be generated by Local Operations 
and Classical Communication (LOCC), entangled states 
must be generated somewhere, and then they have to 
be distributed among different parties, possibly far away 
from each other {transmission) HQ. Once entanglement 
has been distributed, it can then be used immediately or 
stored, waiting for the time to use it [storage). Sys- 
tems which are physical supports for entangled states, 
unavoidably interact with the environment, both during 
transmission and storage times, and therefore undergo 
noisy processes that damage entanglement. For practi- 
cal purposes it is necessary to quantify the amount of 
entanglement that a system loses during time evolution. 

For a pure state p = bipartite entanglement 

between subsystems A and B is unambiguously defined 
as the entropy of entanglement E(\ifj)(ip\) — S(pa) = 
S(pb), where S(pi) is the von Neumann entropy of one 
of the two reduced states, pa = Tr^/o and ps = T^aP- 
The quantification of entanglement for mixed states is a 
much more complicated, open problem 0, 0|. The diffi- 
culty here is rooted in the fact that a mixed state p may 
be decomposed into an ensemble of pure states in an in- 
finite number of ways, p = J2 i Pi\ipi)(4'i\, with pi > 
and ^2iPi = 1. The arbitrariness of the decomposition 
renders any quantification of mixed-state entanglement 
difficult, since some optimization over all possible decom- 



positions is needed. 

In this Letter, we show that any quantification of en- 
tanglement uniquely based on the density operator p 
may lead to the paradoxical result that the entangle- 
ment of bipartite systems can increase under local op- 
erations. Since local operations cannot create entangle- 
ment, the increase of entanglement must be attributed 
to the manifestation of quantum correlations that were 
already present before the application of the local op- 
eration itself. These quantum correlations are in some 
sense hidden, since the density operator formalism does 
not capture their presence. We propose a quantitative 
definition of Hidden Entanglement (HE) and illustrate 
the usefulness of this concept by several examples. We 
show that the presence of HE allows to recover entangle- 
ment between two noninteracting qubits subject to a low- 
frequency stochastic environment by local pulses (acting 
only on one qubit) . We also discuss the relation between 
HE and the occurrence of entanglement revivals, that is 
an up-to-date open issue. Finally, we propose a quan- 
tum information scheme where the presence of HE can 
be revealed. 

Hidden entanglement.— We deal with an ensemble 
A = {{pi, IV^))} °f systems, that is, we know the statis- 
tical distribution of the bipartite pure states {l^i)}; oc- 
curring with probabilities {pi\, so that p = X^PilV'iXV'il; 
but we do not know the state of any individual system 
in the ensemble. The amount of entanglement present in 
the ensemble A can be properly identified as its average 
entanglement 



E av {A) = Y,PiE{\i>i)&i\)). 



(1) 
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If each system in the ensemble evolves during time t un- 
der LOCC, the maximum amount of entanglement of the 
corresponding density operator p(t) can never overcome 
the initial value E av (A). This statement can be proved 
by the following simple argument. Suppose Charlie pre- 
pares a bipartite system in a (possibly entangled) pure 
state of the ensemble A. Then he sends one half of the 
system to Alice and the other half to Bob through noise- 
less quantum channels. Alice and Bob are linked only by 
a noiseless classical channel. Charlie repeats this oper- 
ation N times. Among these, a certain number Mj of 
times Alice and Bob deal with the state \ipj). If Alice 
(or Bob) receives from Charlie the classical information 
about which state he sent each time, in the limit of large 
N Alice and Bob can distill - by only using LOCC - up to 
Mj E(\ipj)(ipj\) maximally entangles states from the -/In- 
states at their disposal Distillable entanglement 
is in fact the entropy of entanglement for pure bipartite 
states. Therefore, the maximum entanglement that Al- 
ice and Bob can distill per pair, by using the classical 
information received from Charlie, is 



N-^oo IV * — ' A — 



(2) 



where limjv^oo -jf- = Pi, that is just the average entan- 
glement of Eq. ([T]) . 

We define the hidden entanglement of the ensemble 
A = {(pi, as the difference between the average 

entanglement and the entanglement 0, H[ of the state 
p = Pi\ipi) {ipi\ quantified by any convex measurement 
E(p) (reducing to the entropy of entanglement for pure 
states), that is 



E h {A) = E av (A) - E(p) = 



(3) 



Thanks to convexity, Eh is always larger than or equal to 
zero. The meaning of the hidden entanglement of Eq. Q 
can be understood as follows: it measures the entangle- 
ment that cannot be exploited as a resource due to the 
lack of knowledge about which state of the mixture we are 
dealing with. Such entanglement can be recovered (un- 
locked) once this (classical) information is provided [lCf . 
There exist several inequivalent measures of mixed state 
entanglement 0, IH ■ Here we consider the entanglement 
of formation Ef(p), which is an upper bound for any bi- 
partite entanglement measure |13j |. so that E av — Ef(p) 
is a lower bound for the hidden entanglement. Ef(p) 
can be readily computed for two-qubit systems via the 
concurrence C(p) 14 1. 

Random local fields.— We first illustrate the con- 
cept of HE by considering a random, local dynamics 
and demonstrating that, under proper conditions, a com- 
plete recovery of the entanglement Ef(p) of the state 



may occur. A basic property of the average entangle- 
ment is its invariance under local unitary transforma- 
tions acting on any system state IV'i(io)) m the ensemble 
A(t ) = {{Pi, \t/ji(t )})} . In particular this is the case 
of evolution in a random local external field [15| induc- 
ing local random unitaries U a (t) (g) Vp(t), with the op- 
erators U a and V/3 acting on the first and on the sec- 
ond subsystem, respectively, and a, ft random variables. 
The average entanglement is conserved by this dynamics, 
E av (A(t)) = E av (A(to)). On the other side, the entan- 
glement of the mixture is only upper bounded by the 
average entanglement: Ef(p(t)) < E av (A{tj), implying 
that a variable (time dependent) HE may exist. This is 
clearly illustrated by the following simple example. 

Let us consider a two-qubit system AB initially pre- 
pared in the maximally entangled Bell state \(f> + ) (we de- 
note the states of the Bell basis as \<jt k ) = ^(|00) ± |11)) 

and IV^) = ^(|01) ± 1 10) )). The time evolution con- 
sists of local unitaries, but we have no complete infor- 
mation about which local unitary is acting. In par- 
ticular, we suppose that the qubit A undergoes, with 
equal probability, a rotation about the x-axis of its 
Bloch sphere, U x {t) — e"" 1 *"'/ 2 , or a rotation around 
the z-axis, U z (t) = e r ia ^ t/2 , while the qubit B re- 
mains unchanged. Hence, the ensemble A at time t is 
A(t) = {(§, (U x (t) ® l B )|</>+}) , (±, (U z (t) ® 1b)I0 + ))}- 
Since we are dealing with random local unitaries, the av- 
erage entanglement is constant in time, E av (A(t)) = 1. 
On the other hand, the entanglement of the state p{t) 
changes in time. Att= ^, p(t) = i|<?T)(<?!>-| + ±|?/> + )(?/> + | 
is separable, whereas at 2t, U x (2t) = U z (2t) = 1a and 
the initial maximally entangled state is recovered. In the 
interval [?, 2i] the entanglement revives from zero to one 
without the action of any nonlocal quantum operation, 
thus apparently violating the monotonicity axiom. The 
ensemble description tells us that at time t the system 
is always in an entangled state {\4>-) or \ip+)), but the 
lack of knowledge about which local operation the system 
underwent prevents us from distilling any entanglement: 
entanglement is hidden, E h (A(t)) = 1 and Ef(p(t)) = 0. 
At time 2i this lack of knowledge is irrelevant since the 
two possible time evolutions result in the identity oper- 
ation 1a and entanglement is recovered, Eh(A(2t)) = 
and Ef(p(2t)) = 1. We notice that entanglement revivals 
have already been observed for a two-qubit system under 
classical environments [lil 17 1 and they have been in- 
terpreted in terms of correlations in a classical-quantum 
state of environments and qubits 17[, but not in terms 
of HE. 

Stochastic low-frequency noise. A fingerprint of 
the existence of HE is the possibility to completely re- 
cover entanglement of a noisy bipartite system by the ac- 
tion of local pulses. Here we consider a simple system con- 
sisting of two noninteracting qubits affected by stochastic 
low-frequency noise. This simplified model captures es- 
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scntial features of several nanodevices whose dynamics is 
dominated by low- frequency noise [l8j . We suppose the 
two qubits are initially prepared in a Bell state \ifo) and, 
for the sake of simplicity, assume that only one of the two 
qubits, namely qubit A, is affected by phase noise (pure 
dephasing), as described by (K = 1) 



■H A (t) = [-fl A a z + e{t)a z + V(t, t)a x }/2 , 



(4) 



whereas qubit B evolves unitarily under Hamiltonian 
T-Lsit). To start with, we suppose that the stochastic 
process e(t) is sufficiently slow to be considered static 
during the evolution time i, with a value randomly fluc- 
tuating from one quantum evolution (quantum trajec- 
tory [19() to the other: e is thus assumed to be a Gaus- 
sian random variable with zero expectation value and 
standard deviation a. V(t,t) indicates an echo tt pulse 
applied at time t and assumed to be short enough to 
neglect the effect of noise during its application. The 
evolution operator during the pulse is e _,<Tx7r / 2 = —ia x - 
Static noise [l8| produces an effect analogous to inho- 
mogeneous broadening in Nuclear Magnetic Resonance 
(NMR) [2(| ■ The system state is described by the quan- 
tum ensemble A(t) = {p(e)de 1 \ip E (t))}, where \f E {t)} = 
f e -if*-H A (t')dt' ®fe- l So HB{t ' )dt '\Lp ). Note that during 
each trajectory the system acquires a random phase. In 
the density operator description of the system, the in- 
formation about the random phase acquired by the sys- 
tem is disregarded by averaging the evolved pure state 
<p e (t)} with respect to the random variable e: p(t) = 
f dsp(e)\tp e (t))(ip e (t)\. Here p{e) is the Gaussian proba- 
bility density function of e. Following Ref. [2l|, we can 
show that a system prepared in a Bell state evolves in a 
mixture p(t) whose concurrence is 



C(p(t)) 




r (7 t 



\a 2 (t-2i) 2 



< t < i, 
t<t<2t. 



(5) 



The average entanglement of the evolved physical ensem- 
ble is E av (A(t)) = 1: during each trajectory, the sys- 
tem state remains in a pure maximally entangled state. 
Instead, the entanglement of formation Ef(p(t)) is ob- 
tained directly by C{p{t)) [l4[ an d it is plotted in Fig. [T] 
In absence of pulses, Ef decays and almost vanishes 
(Ef(p(t)) ~ 0) at times at ^> 1 (red curve). Differently, 
when a local echo pulse is applied at t = t the entan- 
glement subsequently increases (blue curve), reaching at 
t = 2t its initial maximum value Ef(p(2t)) = E^ nax = 1, 
which equals the average entanglement (dashed line). 
The entanglement of formation can also be evaluated 
when the stochastic process e(t) has a finite correlation 
time t entering the autocorrelation function (e(i)e(O)) = 
0- 2 e -|t|/r^ j n tfiis case, significant entanglement recovery 
is still achievable by a local echo pulse, provided that the 
noise correlation time is sufficiently longer than the time 
when the pulse is applied, r ^>t (dotted lines in Fig. [TJ . 




FIG. 1. (Color online) Entanglement of formation Ef(p(t)) as 
a function of the dimensionless time at. The lower (red) solid 
curve corresponds to the free evolution in the presence of static 
noise, the upper (blue) solid curve is the result of the echo pulse 
applied at time at = 4 (indicated by the arrow), from Eq. 
The dashed line at the top is the system average entanglement 
E a v{A(t)) = 1. Dotted curves are the numerical Ej (p(t) ) evaluated 
for a stochastic process e(t) with correlation time ar and with the 
echo pulse applied. From bottom to top: <tt = 20 (cyan curve), 100 
(orange curve), 200 (green curve), and 500 (purple curve). Perfect 
recovery is obtained in the limit t ft — > oo, corresponding to static 
noise (blue solid curve). 



Therefore, the non-Markovian character of noise is cru- 
cial in order to uncover entanglement by local pulses. 

This example shows the possibility to fully recover the 
entanglement E(p) by a local operation. The physical 
mechanism behind this phenomenon is very simple: the 
local 7r-pulse refocuses the different qubit quantum tra- 
jectories restoring at time 2t the qubit A coherence and 
consequently (qubit B evolves unitarily) causing the en- 
tanglement to reappear, with an efficiency depending on 
the correlation time of the stochastic process. Thus, en- 
tanglement can be recovered without any nonlocal con- 
trol. At first sight this result may appear paradoxical, 
since entanglement is by definition a nonlocal resource. 
In fact, here entanglement is not destroyed during the 
time evolution, as indicated by the average entanglement 
which is maximum at any time, E av (A(t)) = 1. Entan- 
glement is instead hidden, Eh increases due to defocusing 
among the different evolutions of the (maximally entan- 
gled) states of A{t). 

From the examples above we argue that a complete 
recovery of entanglement is possible when various mem- 
bers of the physical ensemble evolve differently from each 
other but unitarily. In this case, even though the evo- 
lution of the ensemble averaged density matrix is not 
unitary, no qubit-environment entanglement is generated 
and environment can be described as a classical sys 
tern 



171 ]. Using a terminology borrowed from NMR [22 1 



we speak of incoherent errors, whereas decoherent errors 
arise when the evolution is non-unitary even for a single 
member of the ensemble. For this latter case the average 
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entanglement decays. 

We point out that HE depends on the quantum ensem- 
ble physically giving the system state. The key point is 
that the physical dynamics constrains the evolved den- 
sity matrix of the system to a specific decomposition. In 
the above examples, the physical decomposition of the 
evolved density matrix is always an ensemble of maxi- 
mally entangled states. Had we considered a dynamics 
such as to give, at t = i , a mixture of separable states 
as decomposition of the system density matrix, no local 
operation would have been capable to recover entangle- 
ment (Eh(t) = 0). The dynamics of the system after the 
application of the pulse proves that, at time i, the two 
decompositions arc not equivalent. 

Finally, we remark that the results of the previous ex- 
amples do not violate the monotonicity axiom 0-0]: en- 
tanglement cannot increase under LOCC. The key point 
is that this axiom is fulfilled by all entanglement measures 
provided we consider local operations which are Com- 
pletely Positive, Trace preserving (CPT) maps. Con- 
versely, in the previous examples entanglement recovery, 
from time i to time 2t, is induced by purely local op- 
erations which arc not LOCC, since the corresponding 
density matrix evolution cannot be described by a CPT 
map. To prove this point, it is enough to observe that 
the density matrix p(2t) = p(0), and therefore driving 
the system from the state p(t) to p(2t) is equivalent to 
driving it to p(0)- This operation cannot be described 
by a CPT map since the (CPT) evolution from time to 
time t is not invertible. 

Discussion.— Now we want to show as the phe- 
nomenon of entanglement revivals that we have examined 
is conceptually different from the revivals that a system 
can exhibit due to the interaction with a quantum envi- 
ronment. Let us consider a very simple example, a two- 
qubit system A-B where A resonantly interacts with a 
quantum harmonic oscillator O via a Jaynes-Cummings 
(JC) Hamiltonian 0, 0, while B is isolated. Initially, 
the two qubits are prepared in the Bell state \4>\ B ) and 
O in its ground state \0o)- hi the interaction picture, 
the Hamiltonian is Hao = g(o~+ a + '- a ^)^ where g is 
the coupling constant, cr+ the qubit rising (lower- 

ing) operator, and (a) the oscillator creation (annihi- 
lation) operator. At time i= iz/g the states of A and O 
are swapped with respect the initial state, and the global 

state becomes |CLi} ® \4>bo)' g ivm S Pab = \0a) ® ^b- 
at time t the A-B entanglement is zero, being completely 
transferred to B-O. Any unravelling of the A-B dynam- 
ics gives at this time a quantum ensemble whose average 
entanglement is zero, so that Eh(F) = 0: local operations 
and classical communication cannot recover any entan- 
glement between A and B. Only the subsequent inter- 
action between A and O can gradually restore the A-B 
entanglement thanks to back-action (at time 2t, when a 
new A-O swapping is completed, the initial state is just 
retrieved). Therefore, the entanglement revival is here 



due to the perfect entanglement back-transfer (a nonlo 
cal operation), as well-known in the literature [2ll . l23ll24 | 



This explanation is unsuitable for the examples con- 
sidered above in this paper, where the environment is 
classical and no entanglement transfer is possible. Dur- 
ing the dynamics, the environment acquires only classical 
information about the system A-B and does not entan- 
gle with A or B. Quantum correlations do not leave 
the system A-B but they are simply not accessible due 
to the lack of classical information. Indeed, if at time 
t (Eh(t) = 1) someone provides A-B with the classi- 
cal information about which random unitary the system 
underwent (in the case of random local fields) or about 
which random phase is added to the system states (in 
the case of stochastic pure-dephasing noise) , then all the 
A-B entanglement could be recovered. 

Often, the non monotonous behaviour of the entan- 
glement between independent subsystems is assumed to 
be a signature a non-Markovianity [Hj]. However non- 
Markovian dynamics themselves cannot explain the oc- 
currence of entanglement revivals between independent 
subsystems, since entanglement is a nonlocal resource. 
Non-Markovianity can be rather considered a necessary 
feature for such a dynamical behaviour, as one can see 
from the example we proposed about two qubits affected 
by low frequency noise. It may be fundamental to distin- 
guish the origin of entanglement revivals depending on 
the nature of the environment. Entanglement can be (i) 
transferred forth and back due to some non-local interac- 
tion (quantum environments), or (ii) hidden and recov- 
ered by local operations (classical environments). From 
an operational point of view, the difference is significant: 
in case (i), when entanglement disappears, no local opera- 
tion or classical communication may be helpful to restore 
it; in case (ii) we can act locally to recover entanglement. 

Non-Markovianity has also been recently interpreted 
also as a "flow of information from the environment back 
to the system" [26| . A distinctive feature of the examples 
we have treated (random local fields and stochastic low- 
frequency noise) is that the information is here purely 
classical. 

Quantum information scheme. We now show an 
example implementablc by current technologies in all- 
optical setups [27] , where the presence of hidden entan- 



glement can be detected by sharing private classical infor- 
mation. We assume that Alice prepares Bell's pairs \(fi + ), 
and sends one half of each pair to Bob. Eve disturbs the 
communication with a phase shifter applied half of the 
times, in such a way that Alice and Bob eventually share 
either of the two states \<fi + ) and \<fi~) with equal prob- 
abilities p + = jT = i. Thus, the overall density matrix 
of Alice, Bob, and Eve reads 



PABE = \\4> + )^\ 



1 



\<f>-)(<f>-\®\K-)(K 



(6) 

where |k )} arc two orthonormal states for Eve's 
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system. The reduced density matrix of Alice and Bob, 
Pab = §(|00)(00| + |11)(11|), is separable. Nevertheless, 
if Alice and Bob agree to measure their qubits in the 
computational basis {|0}, |1)}, at the end of the transmis- 
sion of each qubit they share a classical secret bit. One 
can prove this statement by calculating the correlations 
among the systems involved in the communication. We 
measure the correlations between two quantum systems 
Qi and Q2 by the quantum mutual information [if, 0, S| 
I(Qi : Q2) = S(p Ql ) + S(pq 2 ) - S( PQlQ2 ). We then 
obtain that Eve knows which state Alice and Bob share 
(I(AB : E) = 1). However, Alice and Bob share a bit 
of classical information (I(A : B) = 1) and Eve has no 
knowledge about such bit (I (A : E) = I(B : E) = 0). 
This secret communication is possible due to the pres- 
ence of hidden entanglement. In fact, Alice and Bob 
share entangled pairs with a statistical distribution de- 
scribed by the ensemble A = {(1/2, 10+)), (1/2, \<jr))}, 
having E av {A) = Eh{A) = 1. Notice that, if Alice and 
Bob share separable two-qubit states taken from some 
ensemble B, we have Eh (£>) = and it is straightforward 
to show that, each time, Eve knows which state they are 
sharing: private classical communication is not possible 
in this case. 

Conclusions.— We have introduced the concept of 
hidden entanglement based on the quantum trajectory 
description of the system dynamics, and consequently on 
the quantum ensemble description of the system state. 
We have shown that a recovery of entanglement may be 
performed by only local operations, provided that hidden 
entanglement (HE) is present. There is no violation of 
the entanglement monotonicity axiom because no entan- 
glement is destroyed and created in this case: a nonzero 
HE signals a loss of entanglement that is not due to the 
establishment of quantum correlations with the environ- 
ment. Quantum correlations remain within the system, 
but they are not exploitable due to the lack classical in- 
formation. 

We have shown that HE can be of practical relevance 
in solid-state devices which prevalently suffer from low- 
frequency noise, indicating the amount of entanglement 
that can be recovered by standard local pulses, indepen- 
dently of how far apart the qubits are. This allows us, 
differently from other proposals [29j . to avoid resorting 
to non-local gates which may be a much more demanding 
task. The concept of HE can be also applied in the case 
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M HI. In these 



of environment monitoring strategies 
cases, however, the environment cannot be described as 
a classical system and therefore information cannot be 
retrieved by only local operations: "environment mea- 
surements" are also needed to unlock HE. We have finally 
proposed a quantum information scheme where the pres- 
ence of HE can be detected by sharing of private classical 
information. 
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